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Abstract. We investigate the possible bounds which could be placed on alternative 
theories of gravity using gravitational wave detection from inspiralling compact binaries 
with the proposed LISA space interferometer. Specifically, we estimate lower bounds 
on the coupling parameter u> of scalar-tensor theories of the Brans-Dicke type and on 
the Compton wavelength of the graviton X g in hypothetical massive graviton theories. 
In these theories, modifications of the gravitational radiation damping formulae or of 
the propagation of the waves translate into a change in the phase evolution of the 
observed gravitational waveform. We obtain the bounds through the technique of 
matched filtering, employing the LISA Sensitivity Curve Generator (SCG), available 
online. For a neutron star inspiralling into a 10 3 M Q black hole in the Virgo Cluster, 
in a two-year integration, we find a lower bound co > 3 x 10 5 . For lower- mass black 
holes, the bound could be as large as 2 x 10 6 . The bound is independent of LISA 
arm length, but is inversely proportional to the LISA position noise error. Lower 
bounds on the graviton Compton wavelength ranging from 10 15 km to 5 x 10 16 km 
can be obtained from one-year observations of massive binary black hole inspirals 
at cosmological distances (3 Gpc), for masses ranging from 10 4 to 10 7 M Q . For the 
highest-mass systems (10 7 M Q ), the bound is proportional to (LISA arm length) 1 / 2 and 
to (LISA acceleration noise) ^Z 2 . For the others, the bound is independent of these 
parameters because of the dominance of white-dwarf confusion noise in the relevant 
part of the frequency spectrum. These bounds improve and extend earlier work which 
used analytic formulae for the noise curves. 
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1. Introduction and summary 

The Laser Interferometer Space Antenna (LISA) is a gravitational- wave detector being 
designed for launch in the 2010-2015 time frame pQ. Consisting of a triangular array of 
three satellites orbiting the Sun on an Earth-like orbit, it will use laser interferometry 
to open up the low-frequency gravitational- wave window, to complement the high- 
frequency window currently being explored by ground-based interferometers. It is 
expected to be able to observe waves from known binary star systems, from a background 
of white-dwarf binaries, from inspirals of black holes and other compact bodies into 
massive black holes, and possibly from phase transitions in the early universe. 

LISA may also provide new and interesting tests of fundamental physics. In 
previous papers [21 El E] , we showed how observations of waves from inspiralling compact 
binaries could place bounds on alternative theories of gravity, such as theories of the 
scalar-tensor type (e.g. Brans-Dicke theory), or theories with a massive graviton. 

In scalar-tensor theories, the phenomenon of dipole gravitational radiation modifies 
the damping of the binary orbit, and thereby alters the evolution of the phasing of the 
received wave, compared to what general relativity would predict. We showed, for 
example, that, for the inspiral of a neutron star into a black hole of mass 10 3 M Q at a 
distance of 50 Mpc, the lower bound on the coupling parameter uj would be 2.4 x 10 5 , 
with one year of integration prior to the innermost stable orbit [3]. The bound falls off 
with increasing mass. 

In massive graviton theories, the wavelength-dependent propagation speed of the 
waves and the resulting arrival-time offsets also modify the evolution of phasing of the 
received wave, compared to general relativity. For the inspiral of two 10 6 M black 
holes at 3 Gpc, we showed that the lower bound on the graviton Compton wavelength 
X g would be 5.4 x 10 16 km, over four orders of magnitude larger than the bound inferred 
from solar-system dynamics. 

Several developments have motivated us to revisit these bounds. One is the 
dramatic improvement in the solar-system bound on scalar-tensor gravity via tracking 
of the Cassini spacecraft [5|. The new bound is uo > 4 x 10 4 . Another is the 
increasing interest in massive gravity theories, from the viewpoint both of conventional 
modifications of general relativity |E1IZ|, and of multidimensional or brane- world theories 

m 

The third factor is the availability of an improved noise curve for the LISA detector. 
Earlier work employed an analytic approximation to LISA's instrumental noise based on 
the LISA Phase A study [011101, augmented by an estimate of white-dwarf confusion noise 
in the low- frequency regime ^Tj. The improved curves, based on work by Larson et al. 
and Armstrong et al. [121 HB] , additionally take into account the finite propagation time 
of the laser signals during the passage of the gravitational wave, which is embodied in 
a transfer function that modifies the response of LISA at high frequencies. In addition, 
the curves are available online in a "Sensitivity Curve Generator" (SCG) ^3], that 
permits the user to modify the parameters of LISA (arm length, position noise budget, 
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Figure 1. LISA root spectral noise density vs. frequency, including white-dwarf 
confusion noise. Shown is the nominal LISA curve from the Sensitivity Curve 
Generator, and the analytic curve used in i] 

etc) so as to explore the capabilities of different hypothetical LISA instruments. The 
main difference between the analytic curves used earlier and those from the SCG, apart 
from the characteristic oscillations at high frequency arising from the transfer function, 
is that the amplitude noise level of the latter is roughly a factor of 3 higher (see figure 
1). This will have an effect on some of the bounds we report. 

We use the method of parameter estimation via matched filtering to estimate 
bounds on alternative theories. The idea of matched filtering is to cross-correlate a 
theoretical template gravitational waveform against the output of the detector. The 
template may depend on a number of parameters, such as the masses of the stars, and 
parameters associated with the theory of gravity. One can then calculate, for a given 
set of parameters and for a given spectrum of noise in the detector, how accurately the 
parameters inherent in the true signal can be estimated [EH EH H3 CHI UE] • 

For testing scalar-tensor theories, the best system is binary inspiral of a neutron 
star into a black hole. We consider quasi-circular inspiral of non-spinning bodies ending 
at the innermost stable circular orbit (ISCO). Assuming detection with a signal-to-noise 
ratio (SNR) of 10, we find lower bounds on uj shown in figure 2, for various integration 
times. Also shown are representative distances of the sources observed, along with the 
solar-system bound from Cassini. The sharp decrease in the bounds for BH masses 
below a few hundred solar masses results from the fact that the signal falls off the high- 
frequency end of the LISA sensitivity curve before the system reaches the ISCO. The 
bounds shown in figure 2 are only slightly changed from our earlier bounds (figure 1 of 
jlj). For a fixed SNR, the bounds depend on the shape of the noise curve, not on its 
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overall scale (the oscillations from the transfer function tend partially to average out). 
The change occurs in the distance at which a given source can be detected, or in the 
bound obtained from a source at a given distance. Both have decreased by about a 
factor of three, consistent with the shift observed in figure 1. We also find that, for a 
source at a given distance, the bound is relatively insensitive to LISA arm length but is 
inversely proportional to LISA position error. 

For testing massive graviton theories, the best systems are massive black hole 
binaries. Again we consider quasi-circular inspiral of non-spinning bodies ending at 
the innermost circular orbit. For sources at a redshift Z — 1/2, we find lower bounds on 
X g shown in figure 3. The bounds are plotted for various total masses as a function of 
the reduced mass parameter t] = m 1 m 2 /(m 1 + m 2 ) 2 , which varies from zero to 1/4. The 
curves terminate at the small mass ratio end where the SNR drops below 10. For sources 
at a given distance, the bounds on X g are weaker than those reported earlier [3 j by a 
factor of about \/3 (square root because of the quadratic dependence of the effect on 
X g ), again reflecting the higher noise level generated by the SCG compared to our earlier 
formula. For the highest mass systems, e.g. 10 7 x 1O 7 M , the bound is proportional to 
(LISA arm length) 1//2 and to (LISA acceleration noise) -1 / 2 . This dependence becomes 
progressively weaker with decreasing mass, so that for 10 4 x 1O 4 M systems, the bound 
is independent of these parameters. This is because the signals from high mass systems 
reside at the low frequency end of the LISA noise curve, where arm length or acceleration 
noise affect the noise in the expected fashion, whereas the signals from lower mass 
systems reside in the regime where the noise is dominated by white-dwarf confusion 
noise. 
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Figure 3. Bound on graviton Compton wavelength X g vs. 77 = mim2/(mi + t^) 2 , 
for massive black hole binaries at Z = 1/2. 



The bounds shown in figure 3 should be compared to bounds inferred from the 
validity of static Newtonian gravity over large distances, of 2.8 x 10 12 km from solar- 
system dynamics (201 , an d 6 x 10 19 km from cluster dynamics 12 1 j . 

The remainder of this paper provides details. In section 2 we review the method 
of parameter estimation using matched filtering within general relativity. Section 3 
treats scalar-tensor bounds and massive graviton bounds. Section 4 presents concluding 
remarks. 

2. Parameter estimation using matched filtering 

2.1. Estimation in general relativity 

In this section we review briefly the technique used for estimating parameters of 
inspiralling binaries using matched filtering, in the context of general relativity. Further 
details may be found in [TBJ H3 EH d] . In matched filtering, a template consisting of 
a theoretical gravitational waveform is compared to the detector output. The template 
that is a match with an actual signal present with the noise will show a strong correlation 
with the detector output and thereby will "filter" a signal out of background noise. The 
template is represented by h(t), related to the spatial components of the radiative metric 
perturbation far from the source. We adopt the so-called "restricted post-Newtonian 
approximation" for an inspiral orbit that is quasi-circular (that is, circular apart from 
the adiabatic decrease in separation), wherein we approximate h(t) ~ hofye 1 ®^, where 
ho(t) is a slowly varying wave amplitude which depends on the wave polarization, the 
location of the source on the sky, the detector orientation, and the distance to the source. 
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The wave phase, <&(£), is a function of the evolving orbital frequency. Because the process 
of matched-filtering is sensitive mostly to the phasing of the wave, the amplitude h (t) 
is assumed to be given by the lowest-order, quadrupole approximation, while $(£) is 
taken to a suitably high post-Newtonian order. The Fourier transform of h(t) using the 
stationary phase approximation is given by 

f Af- 7 / % e m ^ < f < f 
Hf) = \ ' U < /</max ' (1) 

I, " j J /max ; 

where / max is the largest frequency for which the wave can be described by the restricted 
post-Newtonian approximation, often taken to correspond to waves emitted at the 
innermost stable orbit (ISCO) before the bodies plunge toward each other and merge. 
A useful approximation for this frequency is (we use units in which G — c — 1) 

/max = (e^vrm)- 1 , (2) 

where m is the total mass of the system. After averaging over all angles, the amplitude 
A is given by 

1 A4 5 / 6 

^2/3 D L ' {d) 



where Dl is the luminosity distance of the source, and M. = rf^m is the "chirp" mass. 
In general relativity, the phasing function ^f(f) is given, through 1.5 post-Newtonian 
(PN) order, for spinless bodies, by 

* (/) = 27T/f c - & ^ 



4 

3 

128 



+ u-^ 



1 20 /743 ll \ 2/5 2/3 _ 3/5 ' 

H 1 77 77" 2/ V /3 - l6?T77 _3/5 n 

9 V 336 4 7 ' 



(4) 



where 77 = irAif (77 ~ t; 3 ) and C is formally defined as the phase of the wave at the 
time of coalescence, t c . (Terms through 3.5PN order are known [22] but will not be used 
here.) 

By maximizing the correlation between a template waveform that depends on a 
set of parameters 8 a (for example, the chirp mass Ai) and a measured signal, matched 
filtering provides a natural way to estimate the parameters of the signal and their errors. 
With a given detector noise spectral density S n (f) one defines the inner product between 
two signals hi(t) and /12 (t) by 

(h 1 \h 2 )=2 — - — df, (5) 

Jo b n {J) 

where h\(f) and h^f) are the Fourier transforms of the respective gravitational 
waveforms h(t). The signal-to-noise ratio (SNR) for a given h is given by 

p[h] = {h\h)W . (6) 

One then defines the "Fisher information matrix" T a b with components given by 



dh 


dh\ 
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An estimate of the rms error, A9 a , in measuring the parameter 9 a can then be calculated, 
in the limit of large SNR, by taking the square root of the diagonal elements of the inverse 
of the Fisher matrix, 

A6 a = Vz™, E = r- 1 . (8) 

The correlation coefficients between two parameters 9 a and 9 b are given by 

c ab = £ a6 /V£ aa £ b6 . (9) 

In general relativity, the parameters to be estimated using the above template would 
be (j) c , f t c , \nA4, and In 77, where f is a fiducial frequency characteristic of the detector 
noise spectrum. The method then follows that used, for example, in ^Hj: combining 
equations (JTJ) and (JH) and calculating the partial derivatives dh/d9 a for the four listed 
parameters, we construct the Fisher information matrix using equations © and (J2J). 
We then invert the information matrix and evaluate the errors in the four parameters, 
along with the correlation coefficients, notably between Ai and 77. For the alternative 
Brans-Dicke or massive graviton theories to be considered here, we will add a suitable 
term to the phasing *&(f), dependent upon an additional parameter 9. We will augment 
the dimension of the Fisher matrix by one, and estimate five parameters, along with 
the correlation coefficients between Ai, rj and 9. Since the nominal value of 9 will be 
assumed to be zero (corresponding to uj = 00 or \ g = 00), the error on 9 will translate 
into a lower bound on u> or \ g . 

2.2. Space-based interferometers 

We consider space-based interferometers of the proposed LISA type, with a sensitive 
bandwidth between 10~ 5 and 1 Hz, and an expected noise curve which can be expressed 
in terms of an overall amplitude So, and a function of the ratio x = f/fo'- 

S n (f) = S g(x) . (10) 

In earlier work jSJ E] we adopted an analytic noise curve that included the LISA 
instrumental noise and an estimate of "confusion noise" from a population of galactic 
white-dwarf binaries [§J EH IU| given by 

S = 4.2 x 1(T 41 Hz" 1 , 

f =10- 3 Hz, 

g(x) = VWX~ U/3 + 1 + X 2 /1000 + 313. 5ar (6.398+3.5181og 10 *) _ (11) 

In this paper, we will retain the scaling using So and fo, which is useful for expressing 
analytic results, but for g{x) we will use values obtained from data files available on line 
using the LISA "Sensitivity Curve Generator" (see section l2~3j) . 

In terms of this scaling, the SNR is given, from equations ((TJ), © and (0), by 

p 2 = A\A\ 2 fo 4/3 I(7)/S , (12) 
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where we define the integrals I(q) by 

I(q)=l —-dx, (13) 

where x min = / m i n / fo and x max = f mgx / fo, corresponding to the minimum and maximum 
frequencies over which the detector will integrate . In some calculations, the maximum 
value of /max corresponds to radiation emitted at the ISCO of the system, while in 
others, we can consider the effect of terminating observations sooner than this final 
orbit. The frequency / TO i n corresponds to the gravitational- wave frequency observed a 
time T earlier, where for LISA-type systems, we will assume T < two years. Using the 
quadrupole approximation for radiation damping, which is an adequate approximation 
for this purpose, one can relate the frequencies of gravitational radiation at the beginning 
and end of any time interval T by the expression 

/ 256 T 8/3 \^ 3/8 
Ui = «,(! + — -u f j . (14) 

The luminosity distance Djj to which such sources can be detected can be obtained 
from equations @ and (fH2|) . and related to source masses, detector characteristics and 
the SNR: 

D , 2~.M 5/6 



1/2 




4.81Gpc(4r/) 



2. 3. Sensitivity curves for LISA 

In this paper, we shall adopt sensitivity curves for LISA developed independently by 
Larson et al. and Armstrong et al. |T3]. The two methods are in substantial 

agreement, and the former has been summarized in the Sensitivity Curve Generator 
(SCG), available online [T3|. The sensitivity curves incorporate sources of instrumental 
noise in LISA, such as laser shot noise, acceleration noise, thermal noise, etc., coupled 
with a LISA transfer function which takes into account the effect of the finite time of 
propagation of the laser beams during the passage of the gravitational waves. We assume 
the case where all three arms of LISA are of equal length, and we assume that averages 
over angles and polarizations have been done. The SCG permits various choices to be 
made for LISA instrumental parameters, and has an option to include an estimate for 
confusion noise resulting from a background of galactic white-dwarf binaries jTT] ; this 
background is included in our analyses. 

For the baseline LISA proposal, the noise root spectral density (h per VHz) is 
shown in figure 1, along with the analytic curve used in earlier work. The oscillations at 
the high-frequency end result from the transfer function. Notice that the SCG amplitude 
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Figure 4. LISA noise root spectral density for varying arm length 
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Figure 5. LISA noise root spectral density for varying acceleration noise 



noise level is about a factor of 3 larger than the analytic curve. This has the effect of 
weakening our bounds on uo and X g relative to earlier work. 

In order to study the effect of varying LISA parameters on our tests of alternative 
theories, we vary the arm length, the acceleration noise budget and the position noise 
budget. The resulting curves (including the WD confusion noise) are shown in figures 
4-6. Varying arm length and acceleration noise affects primarily the low-frequency 
noise, while varying the position noise budget affects primarily the high-frequency end. 
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Frequency (Hz) 

Figure 6. LISA noise root spectral density for varying position noise 

Varying other parameters in the SCG, such as laser power or wavelength, or telescope 
diameter, has no sizable effect on the noise curve or on the bounds. 

3. Testing alternative theories of gravity 

3.1. Bounding scalar-tensor theories of the Brans-Dicke type 

One of the most striking differences between scalar-tensor theories and general relativity 
is the prediction of dipole gravitational radiation. The source of such radiation is the 
difference between the self-gravitational binding energy per unit mass between the two 
bodies in a binary system, as encoded in a "sensitivity" s (technically, the sensitivity 
of the body's total mass to variations in the locally measured value of the gravitational 
constant). For black holes, s = 0.5, while for neutron stars, s ~ 0.1 — 0.2, depending 
on the equation of state. In principle, dipole radiation can be a strong effect, because 
it depends on two powers of velocity fewer than quadrupole radiation. This additional 
source of energy flux will alter the inspiral orbit of the binary, and thus will modify the 
evolution of the waveform phasing. 

To sufficient accuracy, this can be taken into account in our phasing function ^(f) 
by adding the following term to the general relativistic formula, equation (J3J): 

S*(f) = A M -5/3 f_l^ 2/V2/ 3 ] (16) 

u ; 128 V 84 w / 
where S = si — s 2 (see 0111 for details). In doing so, we are assuming that lu 1, to 
be consistent with the solar-system bound to > 40,000, so that the 0(1/ u) corrections 
to the general relativistic terms in equation (J3J) may be ignored compared to the dipole 
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1000 


3.94 


18.3 


.000222 


.1034 


203772 


.886 


-.994 


-.929 


5000 


3.78 


11.5 


.000528 


.0246 


58020 


.970 


-.998 


-.954 


10000 


5.27 


12.6 


.000739 


.0174 


31062 


.976 


-.997 


-.957 



Table 1. Estimated parameter errors for 1.4 M© NS-MBH systems in Brans-Dicke 
theory: SNR = 10, integration time is one year prior to the ISCO, neutron star 
sensitivity sjvs = 0.2. 

term of equation (JTfij) . Notice that, compared to the leading term "1" in equation (jlj), 
the dipole term is 0{u~ 2 ^) ~ 0(1 /v 2 ). 

Because s = 0.5 for black holes, binary black holes do not emit dipole radiation at 
all. For neutron stars, s is a relatively weak function of mass, and so for binary neutron 
stars with masses near 1.4 M , dipole radiation is suppressed by symmetry. The only 
promising sources, then, are mixed, such as black-hole neutron-star systems (black-hole 
white-dwarf systems were discussed in 01). 

We now carry out the parameter estimation calculation for the parameters <p c , fot c , 
lnA4, Inn, and b = 5S 2 r) 2 / 5 /48u; (b — ► in the GR limit). We consider a neutron star 
of mass 1.4 M in a quasicircular inspiral orbit around a black hole of mass mBH- We 
adopt the value s^s — 0.2. The results of the parameter estimation for various black- 
hole masses are shown in table 1, and the bounds on u as a function of black hole mass 
and for various integration times are shown in figure 2. For black holes less massive than 
a few thousand solar masses, the bounds could exceed the current solar system bound, 
figure 7 shows, for a 10 3 M black hole, the bound that could be achieved for sources at 
various redshifts, and as a function of the LISA position error budget (2 x 10~ n m/v^Hz 
is the nominal value). Since Ab oc 1/p, the bound on u should be inversely proportional 
to the position error budget, as confirmed in figure 7. The bounds can also be shown 
to be relatively insensitive to LISA arm length or acceleration noise. These results are 
consistent with figures 4-6, since the relevant sources are concentrated at the high- 
frequency end of the LISA noise spectrum, where the noise spectral density is affected 
by position errors but not by arm length or acceleration noise. These bounds for sources 
at a set distance, are roughly a factor 3 weaker than the bounds estimated in as 
expected from the overall shift shown in figure 1. 

3.2. Bounding massive graviton theories 

Notwithstanding statements in the literature forbidding theories of gravity with a 
massive graviton [231 121]; such theories have attracted considerable recent interest, 
from the point of view of conventional modifications of classical GR [23] El U\, and 
from considerations of extra dimensions |8,. To the extent that a candidate massive 
graviton theory merges smoothly with GR in the limit that the graviton mass m g 
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Figure 7. Bound on u vs. LISA position error, for sources at various redshifts. 
Shown are values of the SNR for representative points. 



vanishes or its Compton wavelength X g — > oo, one might reasonably expect corrections 
to GR results for the intrinsic behaviour of a system to be of order (R/X g ) 2 , where 
R is the characteristic size of the system. Since solar system measurements already 
constrain X g > 10 12 km, such corrections can be ignored for inspiralling compact 
binaries. The dominant effect of a massive graviton is via the propagation of the waves: 
lower frequency waves propagate more slowly than higher frequency waves. For binary 
inspiral seen at cosmological distances, this wavelength dependent speed difference can 
result in a significant cumulative distortion of the apparent phasing of the observed 
chirp gravitational-wave signal. The effect is to add the following term to the general 
relativistic formula, equation (jlj): 

3 _ 5/3 / 128 ti 2 DM 2/3 \ 

5V (f) = u b/6 —77 -u 2 ' 6 , (17) 

KJJ 128 V 3 x K l + z ) 1 

where Z is the redshift of the source, and D = (1+Z) f£[a(t) / a(to)]dt, where a(t) is the 
cosmological scale factor, and t e and t are the times of emission and reception of the 
gravitational- wave (see [3] for details). The "distance" D arises from the wavelength- 
dependent propagation of the gravitational wave signal. It is related to the normal 
luminosity distance D L by D/D L = [1 + (2 + Z)(l + Z + VTT~Z)}/5(1 + Zf in a 
spatially flat, matter dominated universe. 

We now carry out the parameter estimation calculation for the parameters <p c , fot c , 
ln.M, ln?7, and f3 = ir 2 DM/ X 2 g (l + Z) ((3 — > in the GR limit). We consider massive 
binary black hole systems without spin, in quasicircular inspiral orbits, for year-long 
integration times leading to the ISCO. The sources are assumed to be at 3 Gpc, with 
a Hubble constant of 70 km/s/Mpc, and a spatially flat universe. The results of the 
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10 7 
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1023 


31.3 


.119 


.0316 


.945 


4.8 x 10 16 


-.979 


-.991 


.997 


10 6 


10 6 


984 


1.65 


.048 


.0050 


.276 


3.1 x 10 16 


-.965 


-.988 


.993 


10 5 


10 5 


871 


.202 


.034 


.0015 


.140 


1.5 x 10 16 


-.953 


-.986 


.987 


10 4 


10 4 


128 


.756 


.286 


.0013 


.441 


0.4 x 10 16 


-.957 


-.988 


.989 



Table 2. Estimated parameter errors for binary MBH systems at 3 Gpc in massive 
graviton theories: integration time is one year prior to the ISCO. 



parameter estimation for various black-hole masses are shown in table 2. The bounds 
on X g as a function of black hole mass for sources at Z — 1/2 are shown in figure 3. 

It is useful to note that, in the large p limit, all errors such as A/3 are inversely 
proportional to the SNR p. Defining B 1 / 2 = ^(S^) 1 / 2 , viewing A/5 as an upper bound 
on (3, and combining this definition with equations (jHJ), © and (|T2j) we obtain an 
expression for the lower bound on \ g : 

Xa> [T5^) [(1 + Z)D L ) f ^ B i,A ■ (18) 

Note that the bound on \ g depends only weakly on distance, via the Z dependence of 
the factor [-0/(1 + Z)D L }^ 2 , which varies from unity at Z = to 0.68 at Z = 1/2. This 
is because, while the signal strength and hence the accuracy fall with distance, the size 
of the arrival-time effect increases with distance. Otherwise, the bound depends only 
on the chirp mass and on detector noise characteristics. 

The resulting bounds are in the range between 10 15 and several xlO 16 km. Figures 
8 and 9 show the dependence of the bounds on LISA arm length L and acceleration 
noise A. For the highest mass systems, the bounds vary as L 1//2 , and A^ 1 ^ 2 because 
the noise varies as L~ l and A, while for the lower mass systems, there is little effect 
because the relevant signal is in the regime where WD confusion noise dominates over 
instrumental noise. 



4. Discussion 



We have used up-to-date LISA sensitivity curves from the SCG to refine the bounds 
that could be placed on alternative theories of gravity. 

For scalar-tensor theory, the bound on uj could exceed that from the best current 
solar-system measurement, using the Cassini spacecraft, by factors of 10 and higher. For 
neutron star inspiral into a 400 M & black hole at 10 Mpc and a two-year integration, 
the bound could reach 2 million. 

On the other hand, a number of alternative bounds on scalar-tensor gravity may 
reach levels competitive with these bounds on a comparable timescale. The gyroscope 
experiment Gravity Probe-B, set for launch in April 2004, anticipates placing a bound 
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LISA Arm Length (10 6 km) 

Figure 8. Bound on graviton compton wavelength X g vs. LISA arm length for sources 
at Z = 1/2; 5 x 10 6 km is the baseline. 




Acceleration noise (10~ 15 m s" 2 Hz" 1/2 ) 



Figure 9. Bound on graviton compton wavelength X g vs. LISA acceleration noise 
for sources at Z = 1/2; 3 x 10~ 15 ms~ 2 Hz -1 / 2 is the baseline. 
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uj > 10 5 via measurement of the geodetic precession of gyroscopes orbiting in the 
curved spacetime around the Earth. Future space optical interferometery missions, 
such as GAIA, could reach comparable levels by measuring the deflection of light to 
microarcsecond precision. GAIA is planned for launch around the same time (2010) as 
LISA. The best hope for a dramatic improvement in u bounds comes from the recently 
analysed binary pulsar system PSR J1141-6545, in which the companion is most likely 
a white dwarf [2S]. Because of the asymmetry between the neutron star and white 
dwarf sensitivities (0.2 vs. 10~ 4 ), dipole gravitational radiation is significant, and a 
measurement of the rate of orbital decay P& in agreement with GR at the one percent 
level could bound ui by as high as 10 6 |23EB]- Such a result could possibly be reached 
in a decade, the same timescale as LISA. 

A major uncertainty in our proposed bound using LISA is the likelihood of 
observation of relatively nearby inspirals of a neutron star into intermediate mass black 
holes. Miller has estimated the rate of inspirals of intermediate-mass binary black 
hole systems in globular clusters. For a 10 M black hole inspiraling into a 100 M 
black hole, Miller estimates yields a rate of one every 250 years, for one-year integrations. 
However one of us [33J has used the SCG to estimate a rate 7 x 10 5 times lower. Whether 
these estimates apply to the neutron-star inspirals needed for scalar-tensor bounds is 
an open question at present. Although the rate may be low from the point of view of 
gravitational- wave astronomy, it is useful to point out that, to test alternative theories, 
a single serendipitous discovery is all it takes: witness the Hulse- Taylor binary pulsar. 

Other methods have been suggested for bounding the graviton mass. We earlier 
showed that observations of binary inspiral using ground-based detectors of the advanced 
LIGO type could place a bound of several times 10 12 km, comparable to the solar-system 
bound [3]. Sutton and Finn showed that, in a simple class of linearized massive 

graviton models, a bound comparable to or slightly better than the solar-system bound 
could be placed on \ g using binary pulsar data. Cutler, Hiscock and Larson [32J studied 
the effect of massive gravity on the observed gravitational-wave phase of known binary 
stars, compared to the orbital phase inferred from light measurements and showed that 
LISA observations could place bounds on \ g from 5 to 50 times larger than the current 
solar-system bound. 
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